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x(t)  -  Ax(t)  +  Bu(t) 


(1.1) 


and  its  discrete  counter  part 


Xj^jft)  *  Ax^t)  +  Bu(t) 


(1.2) 


Uhere  A  and  B  are  constant  matrices  of  appropriate  dimensions  and  x  and 
us  are  time  dependent  vectors,  arise  in  a  wide  variety  of  practical  situations. 
These  Include  design  and  analyses  of  industrial  complexes  such  as  chemical 
plants,  electro-mechanical  machines,  like  motor  cars,  aircraft,  spaceships, 
economics  structures  of  countries,  etc.  Associated  with  these  sy terns  are  many 
interesting  Linear  Algebra  problems.  The  very  well-knowns  are 
STABILITY  AND  INERTIA  PROBLEMS 
CONTROLLABILITY  AND  OBSERVABILITY  PROBLEMS 
POLE  ASSIGNMENT  PROBLEMS 

MATRIX  EQUATIONS  PROBLEMS  (SYLVESTER,  LYAPUNOV,  RICCATI,  etc.) 

At 

In  addition,  there  are  many  other  subproblems  such  as  (1)  evaluating  e 
(2)  determining  Relative  Primeness  of  two  polynolals  or  matrices  (3)  finding 
the  Cauchy  index  of  rational  functions,  etc. 

Computer  solutions  of  these  problems  are  of  utmost  Importance.  Unfortunately, 
only  very  few  viable  sequential  methods  are  available  for  solutions  of  these 
problems.  Moreover,  FAST  (0(n)  or  OCnlogjn)  steps)  PARALLEL  algorithms  for 
these  problems  do  not  exist  at  all  (as  far  as  these  investigators  are  aware). 

The  major  objectives  of  this  project  are  to  develop  fast  computational 
algorithms,  both  sequential  and  parallel,  for  the  problems  mentioned  above 
and,  to  study  parallel  arithmetic  omplexltles  of  these  problems,  i.e.  how  far 
these  problems  can  be  solved  assuming  that  sufficient  number  of  processors  are 
available  [8,  433. 


I 


Up  to  date,  the  following  accomplishments  have  been  made: 


(A)  We  have  developed  FAST  parallel  algorithms  (0(n  log  n)steps  -  0(n  ) 
processors  algorithms)  for  (i)  determinining  controllability  of  the  pair  (A, 8) 
where  A  is  a  companion  matrix  and  B  is  full  (ii)  single-input  pole  assign¬ 
ment  problem  (iii)  stability  problems  for  polynomials  and  (iv)  problems  for 
determining  relative  primeness  of  two  polynomials  and  the  Cauchy  index  of 

a  rational  function.  These  algorithms  have  parallel  efficiency  of  0(r— ^ —  ) 

log  n 

in  most  cases. 

A  desirable  feature  of  these  algorithms  is  that  they  use  only  linear 
algebraic  operations  for  which  parallel  algorithms  already  do  exist. 

(B)  We  have  also  been  able  to  show  (sometimes  reconstructing  the  proposed 
0(n  log  n)  steps  -  0(n  )  processors  algorithms,  sometimes  using  entirely 
different  technqiues)  that  parallel  arithmetic  complexities  of  all  the  control 
problems  mentioned  above  are  upper  bounded  by  0(log  n)  parallel  steps,. 


(C)  We  have  developed  a  FAST  sequential  algorithm  (two  to  three  times  faster 

then  the  best  known  methods)  for  solving  the  symmetric  positive  semidefinlte 

T  T 

Lyapunov  matrix  equation  XA  +  A  X  «  C  C . 

(D)  A  FAST  direct  sequential  method  for  finding  the  eigenvalue  distribution 
of  a  matrix  inside  and  outside  various  regions  of  the  complex  plane  including 
half  planes,  shifted  half  planes,  hyperbolas,  sectors,  quadrants,  imaginary 
axis,  regions  contained  within  two  straight  lines  that  pass  through  the  origin 


(E)  We  have  just  developed  a  parallel  algorithm  for  computing  the  zeros  of 
a  polynomial  whose  zeros  are  all  real  and  distinct.  This  algorithm  requires 
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